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Abstract 

A statistical analysis of the solar neutrino data is presented assuming the solar 
neutrino deficit to be resolved by the resonant interaction of the neutrino magnetic 
moment with the solar magnetic field. Four field profiles are investigated, all exhibiting 
a rapid increase across the bottom of the convective zone, one of them closely following 
the requirements from recent solar physics investigations. First a ’rates only’ analysis 
is performed whose best hts appear to be remarkably better than all fits from oscilla¬ 
tions. A global analysis then follows with the corresponding best hts of a comparable 
quality to the LMA one. Despite the fact that the resonant spin havour precession does 
not predict any day/night effect, the separate SuperKamiokande day and night data 
are included in the analysis in order to allow for a direct comparison with oscillation 
scenarios. Remarkably enough, the best ht for rates and global analysis which is com¬ 
patible with most astrophysical bounds on the neutrino magnetic moment is obtained 
from the prohle which most closely follows solar physics requirements. Allowing for 
a peak held value of 3 x 10®G, it is found in this case that Am^i = 1.45 x 10“®eF^, 
Pi, = 3.2 X (65%CL). The new forthcoming experiments on solar neutrino 

physics (Kamland and Borexino) will be essential to ascertain whether this fact is 
incidental or essential. 
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1 Introduction 


Whereas the apparent anticorrelation of the neutrino event rate with sunspot activity 
claimed long ago by the Homestake collaboration [Q remained unconhrmed by other ex¬ 
periments i, i, i and theoretical analyses [^, the magnetic moment solution to the 
solar neutrino problem is at present an important possibility to be explored in the quest for 
an explanation of the solar neutrino deficit. This is the idea, originally proposed by Cisneros 
P], and later revived by Voloshin, Vysotsky and Okun [0, that a large magnetic moment 
of the neutrino may interact with the magnetic field of the sun, converting weakly active to 
sterile neutrinos. It now appears in fact that this dehcit is energy dependent, in the sense 
that neutrinos of different energies are suppressed differently. In order to provide an energy 
dependent deficit, the conversion mechanism from active to nonactive neutrinos must be 
resonant, with the location of the critical density being determined by the neutrino energy. 
Thus was developed the idea of the resonance spin flavour precession (RSFP) proposed in 
1988 [Q. It involves the simultaneous flip of both chirality and flavour consisting basically 
in the assumption that the neutrino conversion due to magnetic moment and magnetic field 
takes place through a resonance inside matter in much the same way as matter oscillations 
0 . A sunspot activity related event rate in a particular experiment would hence imply that 
most of the neutrinos with energies relevant to that experiment have their resonances in the 
sunspot range. However, the depth of sunspots is unknown (they may not extend deeper 
than a few hundred kilometers) and the observed held intensity is too small in sunspots to 
allow for a signihcant conversion. The anticorrelation argument has therefore lost its appeal 
for several years now. 


Despite the absence of the anticorrelation argument, several main reasons may be in¬ 
voked to motivate RSFP and investigating its consequences for solar neutrinos. In fact both 
RSFP and all oscillation scenarios indicate a drop in the survival probability from the low 
energy (pp) to the intermediate energy neutrino sector (^Be, CNO, pep) and a subsequent 
moderate rise as the energy increases further into the ®B sector. The magnetic held prohles 
providing good hts to the event rates from solar neutrino experiments typically show the 
characteristic of a sharp rise in intensity at some point in the solar interior, followed by a 
progressive moderate decrease . This is in opposite correspondence with the energy 


dependence of the probability in the sense that the strongest held intensities correspond 
to the smallest survival probabilities. Hence RSFP ohers a very suggestive explanation for 
the general shape of the probability, which naturally appears as a consequence of the held 
prohle. On the other hand, from solar physics and helioseismology such a sharp rise and 
peak held intensity is expected to occur along the tachoclyne, the region extending from 
the upper radiative zone to the lower convective zone, where the gradient of the angular 
velocity of solar matter is diherent from zero |1^, |]R]. Furthermore, it has become clear 
0 . 0.0 that RSFP provides event rate hts from the solar neutrino experiments that are 
remarkably better than oscillation ones [|T^. Finally, there are recent claims in the 


literature for evidence of a neutrino hux histogram 0 containing two peaks, an indication 
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of variability pointing towards a nonzero magnetic moment of the neutrino. 

This paper is divided as follows: in section 2 a statistical analysis of the four rates 
(Chlorine, Gallium, SuperKamiokande and SNO) on the light of the standard solar model 
(BP 2000) and the most recent data 0, I, I, i 0 is presented. The solar neutrino 
dehcit is assumed to be originated from nonstandard neutrinos endowed with a magnetic 
moment interacting with the solar magnetic field. Four solar held prohles, all obeying the 
general features described above, are analysed. The hrst three were used in a previous 
paper 13 at the time SNO results were non-existant. Special attention should be payed 
to prohle (4), the one most closely following the requirements originated from solar physics 
investigations [|^]. All Gallium measurements (SAGE, Gallex and GNO) are combined in 
one single data point, so that the number of d.o.f. is 2: 4 experiments - 2 parameters. These 
are the mass square difference between neutrino havours Amh and Bo, the value of the held 
at the peak, investigated in the ranges (0 — 10“^)eG^ and (0 — 300)/cG respectively. This 
Am^ range is related to neutrino resonances in the region from the upper radiation zone to 
the solar surface, that is, the whole region where a signihcant magnetic held is expected. 
The upper bound of 300 kG for the peak held {Bq) at the base of the convection zone is 

, [^. The neutrino magnetic moment is taken throughout 
This means that, since the order parameter is the product of by 


suggested by most authors 
at fiu = 

the magnetic held, a suitable to ht to solar neutrino data may not allow a value of Bq 
well above lO^G, as it would imply a value of /ij, well above 3 x in clear conhict 


with astrophysical bounds on ^ 3 . The global analysis performed in section 3 (rates + 
SuperKamiokande spectrum) therefore excludes all hts with Bq > lO^G. The border line 
cases (hts 2.2, 3.1 with Bq = 1.18 x lO^G, Bq = 1.09 x lO^G) were included. In order 
to allow for a direct comparison with the oscillation cases |^, the day/night 

data from the SuperKamiokande collaboration with 1258 days is included in this global 
analysis, despite the fact that no day/night ehect is predicted by RSFP. Further, since the 
information on the SuperKamiokande total rate is already present in the hux of each spectral 
energy bin, this rate does not enter the analysis, thus following the same atitude as other 
authors [^, [l^, |^. The number of d.o.f. is therefore: 3 rates + 19 x 2 spectral points - 2 


parameters = 39. The results are presented in terms of best hts and allowed physical region 
(95 %GL and 99 %GL) contours on the Aml^, Bq plane. A small discussion and conclusion 
section then follows. 
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Rates only: Chlorine, Gallium, SuperKamiokande 
and SNO 


All profiles studied generically satisfy the requirement of a sharp rise near the bottom of the 
convective zone with a subsequent smoother decrease. The hrst three were introduced in 
a previous paper where separate analyses of rates and spectrum were made (no global 
analysis) at a time when SNO results did not exist: 

Profile (1) [equilateral triangle] 


B = 0 
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B = B, 
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where x is the fraction of the solar radius, xr = 0.70, xc = 0.85 and all units are in Gauss. 
Profile (2) 

B = 0 , X < Xr (4) 
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with Xr = 0.65, Xc = 0.75. 
Profile (3) 

B = Bi 
B = 


B = 2.16 X lO'* 
a;-0.75\2' 
0.04 ) 
Bo 


X < 0.7105 
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(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 


cosh 30(x- 0.7483) ’ 

with Bq = 0.998i?i. 

The fourth prohle is chosen so as to more closely agree with solar physics observational 
requirements |]^. In fact a recent detailed study of the solar internal rotation has lead to 
the estimation of a large scale magnetic held. An upper limit of 300 000 G is derived for a 
held near the base of the convection zone, a value already claimed by previous authors |]^ 
at this particular location. Also a held intensity of approximately 20 000 G at a depth of 
30 000 km below the surface is expected. The simplest prohle satisfying such features is 
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Profile (4) 
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( 10 ) 

( 11 ) 

( 12 ) 

(13) 


B = 2.10* +{x- 0.957) 
with Xr = 0.65, Xc = 0.713. 

The possibility of this relatively high intensity (20 000 G) at the small depth of 30 000 
km corresponding to another peak has been known for some time to be inconsistent with 
the data. These seem to be consistent only with a held intensity decreasing monotonically 
from the bottom of the convective zone to the surface 0- 

The ratios of the RSFP to the SSM event rates Rcaci dehned as before [^, the 
SuperKamiokande one is 


TDth _ 

^SK — 


/f “ dE, /“ dE'JiE',, E,) Ji“ dEd>mP(E)d^ 3- (1 - P(J!))^ 


/I™ dE, S,Z dKf(E',, E,) /|" dE4>(E)d^ 


(14) 


Here 4>{E) is the SSM neutrino hux (hep+®8) and the quantity f{E'^,Ef) is the energy 
resolution function |2^ of the detector in terms of the physical {Ef) and the measured {Ef) 
electron energy {E^ = T + rUe). The lower limit of E^ is the detector threshold energy 
(8em = Eeth with 5.0MeV) and the upper limit is E^max = 20Mel7 |^. For the lower 


I0| and upper [|J] integration limits of the neutrino energy, one has respectively 
T' + 


Em 


, Em = IhMeV {i — B) , Em = IS.SMeH {i = hep). (15) 


The weak differential cross sections appearing in equation (17) are given by 


d(7\Y Gf^TTI^ 


27r 


[igv + grf + [gv - gA]* (l “ (^) “ ^9v‘* - gA)—^\ 


(16) 


for z/gC scattering, with gv = \ + 2sin‘^9w, Qa = h- For z/^e and v^e scattering. 


d(7\Y Gp'^Tne 


\ 2 


[{gv - gA)"^ + {gv + 5 'a)^ ~ e) ~ ~ 


27r 


(17) 


with gv = —\ + 2siri^9w, gA = —|- 

The fourth quantity to be evaluated is the ratio of the charged current event rate by the 
standard charged current event rate in the SNO experiment {rcc)- As recently announced 
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by the SNO collaboration rcc = 0.347 ± 0.029 |^. By definition rcc = Rcc/Rcc where 
Rcc is given by 


Rcc = / dTf 

'jTrn 


dEf(E)P(E) 

J Q 


rE-\-nie Q 

/ dE,^{E,E,)R{E,^^^,E,) (18) 

/me dH/p 


and Rcc stands for the same quantity with the replacement P{E) 1. In (18) Eeff denotes 
the measured total electron energy E^ff = + rrif. with a threshold given in terms of the 

kinetic energy, = Q.75MeV. The function R{Ef,^^^, Eg) is the detector energy response 


19| and Q = lAAMeV with the rest of the notation being standard. The CC cross section 


was taken from K. Kubodera’s homepage, 


The partial event rates for each neutrino component in each experiment were taken from 
p!8| and the solar neutrino spectra from Bahcall’s homepage |Q. The contribution of the 
hep flux to both the Gallium and Homestake event rates was neglected. The analysis for 
the ratios of event rates and electron spectrum in SuperKamiokande was done following the 
standard procedure described in [|. The validity of this procedure and alternative ones 

m 


for solar neutrinos is discussed in refs. 


m- 


The ratios of event rates to the SSM event rates, both denoted by Rj in the following, 
were calculated in the parameter ranges = (0 — 10“^)el7^, i?o = (0 — 30) x lO^G for 

all magnetic field profiles and inserted in the definition. 


ii j'2=i 
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(19) 


with (Ga=l, Gl=2, SK=3, SNO=4). The measured rates Rj^^^ in this equation can be read 
directly from table I. Given the four data points and the number of parameters in the fit, the 
number of d.o.f. is 2. Within the above parameter ranges for Am^i and Bq, 19 local minima 
of xlates were found. They are displayed in table II along with the respective values of the 
goodness of fit (g.o.f.). As mentioned in the introduction, a working value of 
was chosen. Since the order parameter is the product of the neutrino magnetic moment by 
the magnetic field which can be as large as 3 x lO^G at the peak, a value Bq = lO^G, for 
example, at the local minimum, means the existence of a solution with 
In order to remain consistent with most astrophysical bounds , it is therefore unadvisable 
to consider those fits for which Bq is well above lO^G. A close inspection of table II shows 
that all profiles give reasonable or good quality fits for a peak field value up to ~ lO^G. 
Based on this criterion and selecting therefore fits 1.1, 1.2, 2.1, 2.2, 3.1, 4.1 (g.o.f. 34, 34, 
38, 25, 42, 51% respectively) it is seen that one obtains comparable fits in RSFP to those 
in the oscillation cases |T5|, pT| , p2| . In fact the best oscillation solution from the ’rates-only’ 
analysis is the VAG one 


1^ for which the g.o.f. is 52%. | Interestingly enough, the best 


tPere only the main definitions and differences are registered. For the calculational details we refer the 


reader to ref.| W 


t There is however a minor difference in criterion, since these authors considered the SAGE and 
Gallex/GNO rates as two separate ones, hence 3 d.o.f. in their analysis (see table 4 of ref. |l^). 
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of all fits investigated here is 4.1 from profile 4, the one in the closest agreement with solar 
physics observational requirements |^. Whether this remarkable fact is coincidental or 
essential is of course too early to ascertain. 

A global analysis of the selected hts 1.1, 1.2, 2.1, 2.2, 3.1 and 4.1 will now be performed 
in section 4. 


3 Global Analysis 


In addition to the rates, the quantity to be investigated now is the SuperKamiokande spec¬ 
trum 




dE, SZ dE'J(E'„ E,) /I" dEUUEmEZ-lf + (1 - P(E)) 


\Y 

~St^ 


E. ssr dEe SZ dE'J(E',, E,) /I" dEME)^ 


( 20 ) 


with j = 1,38. Following other authors O, M ^ and since the information on the 


SuperKamiokande total rate is already included in the flux of each spectral energy bin, this 
rate will now be disregarded. In order to allow for a straight comparison with other solar 
neutrino physics analyses the day/night spectrum is also included although no day/night 
difference is predicted by RSFP. Hence the definition is as follows 


38 

jl J2 = l 


The measured rates (j = 1,..., 38) in this equation can be directly read from table II of 
ref. 0]. In performing the global fit one has 41 experiments (3 rates and 38 spectral data) 
and two free parameters, namely the mass squared difference between neutrino flavours 
and the peak held value Bq, hence 39 d.o.f.. For the global 


X gl 


2 4-2 
X rates ' X spectrum' 


( 22 ) 


The local minima of x^i were investigated for each of the selected four fits from the previous 
section (1.1, 1.2, 2.1, 2.2, 3.1, 4.1). The allowed physical regions and stability of these fits 
were next investigated. The analysis proceeds in terms of the two parameters Arn^i, Bq 
hence two degrees of freedom, and the results are presented in terms of 95 and 99%CL 
contours (Ax^=5.99 and 9.21 relative to the local best fit). Global fits 1.1 and 1.2 of profile 
1, whose contours merge, and 3.1 of profile 3 are shown together in fig. 1. Global fits 
2.1, 2.2 of profile 2, whose contours also merge, are shown in fig.2, while 4.1 of profile 4 
and its contours is shown in fig. 3. In table III the parameter values for each of these 
global fits including the values of the g.o.f. are also shown. It is seen that all 

g.o.f. are of the same order, with fit 4.1 being, interestingly enough, the best, although 
by a narrow margin. Gomparing with the best global fit from oscillations (LMA solution) 
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I5|, |I^, the g.o.f. of fit 4.1 is similar but slightly smaller (64.7% compared to 67% with 
hxed hep flux) or reasonably smaller (compared to 85% with free hep flux). Since the 
hep flux was considered hxed in the present analysis, only the hrst comparison should be 
held as valid. Owing to the small and probably not physically meaningful day/night effect 
(1.29cr) observed in the SuperKamiokande p[, spectrum and since RSFP predicts no such 
effect, it is to be expected that global hts in RSFP could be marginally worse, although of 
comparable quality. The energy spectrum corresponding to global ht 4.1 is shown in £g. 4 
superimposed on the SuperKamiokande data. Allowing for a peak held value Bq = 3x lO^G, 
the respective value of the magnetic moment, as seen from table III, is fii, = 3.2 x 
with Am^i = 1.45 x 10“®eR^. 


4 Concluding Remarks 


The analysis of prospects for the magnetic moment solution to the solar neutrino problem 
reveals both rate and global hts of reasonable or good quality. Previous analyses on RSFP 
BOHl, unambiguously indicated from the part of the data on rates a preference for solar 
held prohles with a steep rise across the bottom of the convection zone in which vicinity 
they reach a maximum, followed by a more moderate decrease up to the solar surface. 
They showed a quality of rate hts alone which is not shared by the oscillation rate hts 
ITSl . Interestingly enough, this class of prohles is the most consistent one with solar physics 
and helioseismology |^, ]K|. A further prohle (4) was introduced in the present paper 


which is the simplest most closely following the present requirements from solar physics 
T3|. It is important to note that, as found in the present paper, this is the prohle giving 


the best of all RSFP hts, both for the rates and global. For this (global) ht one hnds that 
Am|i = 1.45 X and jj,,^ = 3.2 x (g.o.f. 64.7%) if one allows for a solar held 

at the bottom of the convective zone Bo = 300 OOOG. An analysis of spin havour precession, 
both resonant and non-resonant, was also performed in a recent paper |^, whose authors 
limit themselves to a ’rates-only’ analysis with a single prohle. They hnd a better g.o.f. 
for the resonant (39%) than for the non-resonant case (28%), however substantially smaller 
than the best one obtained here with prohle 4. 


The forthcoming neutrino experiments Borexino |29| and Kamland pO will be crucial 


to test the validity of the RSFP hypothesis against oscillations. In fact, at the present stage 
both the LMA and RSFP solutions give global hts of equally good quality and, as far as 
rates only are concerned, RSFP solutions provide the best of all hts. Furthermore, owing 
to the typical shape of the RSFP survival probability [|^ , which is quite diherent from the 
LMA, LOW and VAC ones at the intermediate energies [M|, the RSFP solution predicts 

event rate than these three 
It is quite likely that the 


for Borexino a clearly smaller neutrino-electron scattering [32 


16 


oscillation solutions whose prediction is essentially the same 
Kamland experiment will either bring conhrmation of the LMA scenario, or either disprove 

El 


it, along with the less favoured solutions, LOW and VAC 


Consequently, in case these 
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are excluded by Kamland, the expectation for Borexino will naturally be that of a strong 
reduction relative to the standard event rate, therefore providing evidence for RSFP, unless 
a clear day/night effect is observed. 


In the present analysis only time averaged data were considered and a htting was made to 
a time constant prohle ’buried’ in the solar interior. If, on the contrary, the active neutrino 
flux turns out to be time dependent, a situation most likely to be interpreted through 
the magnetic moment solution with a time dependent interior held, the present approach 
is obviously inadequate. Averaging the event rates over time implies disregarding possible 
information in the data which otherwise is available if different periods of time are considered 
T7|. The robustness of such a procedure will greatly improve with the accumulation of more 


data. 
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Experiment 

Data 

Theory 

Data/Theory 

Reference 

Homestake 

2.56 ±0.16 ±0.15 

7.7±li 

0.332 ±0.05 


i 


Ga 

74.7 ±5.13 

129±^ 

0.58 ±0.06 

l.i 

SuperKamiokande 

2.32 ±0.085 

5.05±J;° 

0.459 ± 0.005±°;°J® 


i 


SNO 

1.75 ±0.15 

5.05±J;? 

0.347 ±0.029 


19 



Table I - Data from the solar neutrino experiments. Units are SNU for Homestake and 
Gallium and for SuperKamiokande and SNO. The result for Gallium is the 

combined one from SAGE and Gallex+GNO. 


Fit 



B, 

(G) 

xLes/2 d.o.f. 

g.o.f. 

1.1 

7.04 

X 

10 - 

-9 

4.2 

X 

10 ^ 

2.15 

34.1% 

1.2 

7.35 

X 

10 - 

-9 

6.5 

X 

10 ^ 

2.17 

33.8% 

1.3 

8.38 

X 

10 - 

-9 

1.35 

X 

10 ® 

1.32 

51.7% 

1.4 

8.62 

X 

10 - 

-9 

1.67 

X 

10 ® 

1.81 

40.4% 

1.5 

1.0 : 

X 

10 - 

8 

2.34 

X 

10 ® 

1.12 

57.1% 

2.1 

1.28 

X 

10 - 

-8 

9.6 

X 

10 ^ 

1.96 

37.6% 

2.2 

1.29 

X 

10 - 

-8 

1.18 

X 

10 ® 

2.75 

25.3% 

2.3 

1.37 

X 

10 - 

-8 

1.72 

X 

10 ® 

1.51 

47.1% 

2.4 

1.38 

X 

10 - 

-8 

1.96 

X 

10 ® 

2.37 

30.6% 

2.5 

1.43 

X 

10 - 

-8 

2.47 

X 

10 ® 

1.12 

57.3% 

2.6 

1.48 

X 

10 - 

-8 

2.73 

X 

10 ® 

2.19 

33.5% 

2.7 

1.59 

X 

10 - 

-8 

3.23 

X 

10 ® 

0.986 

61.1% 

3.1 

1.45 

X 

10 - 

-8 

1.09 

X 

10 ® 

1.73 

42.2% 

3.2 

1.87 

X 

10 - 

-8 

1.83 

X 

10 ® 

2.15 

34.1% 

4.1 

1.43 

X 

10 - 

-8 

9.8 

X 

W~ 

1.34 

51.2% 

4.2 

1.41 

X 

10 - 

-8 

1.31 

X 

10 ® 

3.41 

18.2% 

4.3 

1.49 

X 

10 - 

-8 

2.04 

X 

10 ® 

1.17 

55.8% 

4.4 

1.49 

X 

10 - 

-8 

2.38 

X 

10 ® 

2.74 

25.4% 

4.5 

1.58 

X 

10 - 

-8 

3.05 

X 

10 ® 

1.10 

57.8% 


Table II - Rate fits: local minima of xlates prohles 1, 2, 3, 4 (eqs.(l)-(13)). Only the 
hts 1.1, 1.2 (profile 1), 2.1, 2.2 (prohle 2), 3.1 (prohle 3) and 4.1 (profile 4) are considered 
for the global analysis (rates + day/night spectrum). The remainder correspond to a value 
of Bq largely exceeding lO^G. See the text for details. 
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Fit 

Am2i(eF^) 

Bo{G) 

xli/39d.o.f. g.o.f. 

1.1 

6.80 X 10-y 

4.0 X 10^ 

36.2 

59.7% 

1.2 

7.07 X 10-9 

6.8 X 10^ 

36.3 

59.2% 

2.1 

1.27 X 10-® 

9.4 X 10^ 

36.0 

60.9% 

2.2 

1.26 X 10-® 

1.22 X 10® 

38.4 

49.5% 

3.1 

1.41 X 10-® 

1.04 X 10® 

35.4 

63.6% 

4.1 

1.45 X 10-® 

9.6 X 10^ 

35.1 

64.7% 


Table III - Global fits with the three rates (Homestake, Gallium, SNO) and the Su- 
perKamiokande day/night spectrum. Fit 4.1, besides being the best, although by a small 
margin, is related to a solar held prohle which agrees most closely with the requirements 
from a solar physics analysis (details in the text). 


Figure captions 

Fig. 1. Global hts 1.1, 1.2 (lower left, corresponding to prohle (1)) and 3.1 (upper right, 
corresponding to prohle 3) and their 95% (dashed) and 99% GL (solid) contours. See also 
table III. 

Fig. 2. Same as hg.l for global ht 2.1, 2.2 (prohle 2). See also table III. 

Fig. 3. Same as hg.l for global ht 4.1 (prohle 4). See also table III. 

Fig. 4. The theoretical prediction of the SuperKamiokande recoil electron spectrum (solid 
line) superimposed on the data set (1258 days) for global ht 4.1 (prohle 4). This is also 
the best ht found in the analysis, for which x%obai — g.o.f.=64.7%) for 39 d.o.f.. 
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